﻿ 4 The discrete Fourier transform (DFT) 1 Application goal We study the discrete Fourier transform (DFT) and its applications: spectral analysis and linear operations as convolution and correlation We use Matlab as a simulation environment 2 Definitions and properties of the DFT The DFT is a Fourier transform applied to a input signal that is discrete Unlike the discrete-time Fourier transform (DTFT), it only evaluates enough frequency components to reconstruct the analyzed signal of finite support ]^ Consider a sequence :x→ such that [0xn= for nN We define its ] 12 DFT of the sequence :x→]^ by: 22Nπ jk n− ⎧ N [] ⎪ ∑ xne N k N≤≤ 12, Xk= (1) [] = 1nN⎨ ⎪ 0, otherwise ⎩ where 1NN N= (2) 21−+ The periodic signal [ xn of period N is obtained from [xn extended by periodicity: ]] ∞ [] xnxnmN=− (3) [ ] ∑ m=−∞ Denote: 2π j− New= N (4) Consider the signal: kn (5) [] [] ∑ NXkxnw= nN= means the sum is made for N consecutive values of the signal [xn (one Here nN= ] period) [ Xk is periodic, of period N, and: ] [ ≤≤ (6) ][] XkXk NkN= 12, The original signal is reconstructed from its DFT using: 2N 1 kn − [] (7) [] ∑ Nxn X kw= N= 1kN 1 kn− (8) [] [] ∑ Nxn X kw= kNN= When we use (8), it results that: [ =−−− (9) ][][][] {} 12xnxn nN nNσσ nσ is the discrete-time unit step) ([ ] 1 Consider the sequences [ ][] ixn with the corresponding DFT’s iXk, i=1,2 Main properties of the DFT, [ ][] iXk, for the corresponding periodic signals ixn, are: a) Linearity + ∈ (10) [ ][][][] 12 1 2, a,b Cax n bx n aX k bX k+↔ b) Time shifting 0kn [ ↔∈ (11) ][] 00, Zxn n w X k n− c) Modulation in time 0kn− −∈ (12) [ ][] 00, Zwxn Xkk k↔ d) Circular convolution Circular convolution is denoted with It is convolution for periodic signals, with the same period N, whose result is another sequence, periodic of period N: []xnxmxnm=− (13) [ ] [ ] 12∑ mN= The DFT is: (14) [ ][][][] 12 12xnxn XkXk↔⋅ e) Modulation in time 1 [] (15) [][] [] 12 1 2xnxn Xk Xk⋅↔ N f) Parseval’s theorem 221 xnXk= [] (16) [] ∑∑ nN kNN== 3 Fast Fourier algorithm for DFT computation The expansion of digital signal processing began with the development of the algorithm for the computation of the DFT: the Fast Fourier transform (FFT) This exploits some properties of the complex exponentials The FFT implementation is efficient when the length of the sequence in time or the period is a power of 2 If this requirement is not fulfilled, the interval can be extended until it becomes a power of 2, by adding null samples 4 Applications for DFT 4 1 Spectral analysis of continuous-time signals This is based on the sampling theorem for finite energy signals Consider a continuous- xt, with the length time signal, () ( (17) supp 0,xtT⊂ ) [] {} and the signal resulted by periodicity (see Fig 1): ∞ () (18) ( ) ∑ TxtxtnT=− n=−∞ 2 Fig 1 Analyzed signal (a); Periodized signal (b) Xω the Fourier transform of the signal (xt and assume ()xt can be Denote by () ) T expanded in a Fourier series with a finite number of terms: 2Nπ jkt T () (19) ∑ kTxt ae= kN=− xt with the step: We sample () T T= S (20) 21N+ and consider the discrete-time signal (21) [ () ] sxnxnT= with the DFT X[k] We can see that, 12π ⎛⎞ (22) [] kaXk Xk=+ ⎜⎟ TT⎝⎠ meaning the elements of the DFT are proportional with the samples of the continuous- time Fourier transform of signal ( xt The signal (xt can be reconstructed from its )) samples (21), if the sampling step satisfies (20) In reality it is less likely that the sampling fulfills the above mentioned condition, or it is impossible In this case, the expansion of ( ) Txt contains an infinite number of terms: 2π jkt∞ T () (23) ∑ Tkxt ce= ∞ k=− Thus, aliasing of coefficients appears after sampling: 3 ∞∞ 12π ⎛⎞ (24) ==++ [] () () ∑∑ ⎜⎟ (2 1)21kNpXk c X k N p++ ⎝⎠ ppTT=− ∞ =−∞ >, then the values TX[k] are If the Fourier transform x(ω) is small for 2/NTωπ accurate estimates of the samples of this transform in k2π/T, k=[-N,N] In some spectral analyzers with FFT, linear interpolation of samples X[k] is used If the signal () xt has a large support, we apply the technique of the short time Fourier transform (STFT), see Fig 2 Basically, segments of the signals are analyzed The segments are obtained by multiplication of the signal with a window The simplest case is the rectangular pulse, but in order to avoid Gibbs phenomenon, other signals are used such as Hamming, Gauss, etc Fig 2 (a) Original signal; (b) Rectangular pulse; (c) Analyzed signal 4 2 Computation of convolution Consider two sequences [xn and [xn: ]] 12 > [] xnnNnN= [] xnnNnN=<2340 for and For the computation of convolution of the sequences, the following steps are performed: Obtain signals [xn and [xn by periodicity • ]] 12 1NN N N N≥− −+ (26) 2143+ that is a power of 2; thus the convolution of the two signals can be recovered from the circular convolution of the periodic signals Compute [] [] []XkXkXk= • 12 Compute inverste DFT, [xn, the circular convolution of [xn and [xn • ]]] 12 4 • Compute x[n]: {[] [ ] [ 1]xn n N N n N Nσσ= −− − −−− (27) } 13 2 4 4 3 Computation of correlation Consider the same signals as before (eq (25)), but real-valued Their intercorrelation function is: ∞ =−∗ [] [ ] [ ] [ ] []Rnxmxmnxnxn=+ (28) 1212 1 2XX∑ m=− ∞ To evaluate this function, same steps are made, using convolution betwen the first signal reversed in time, and the second signal 5 Practical part 5 1 Consider the following signals: (a) () 1xt, din fig 3 a (rectangular pulse) (b) () 2xt, din fig 3 b (triangle pulse) σ=−−⎡⎤ (a sine) (c) () ( ) ()() 30sinxtA t t tTωσ ⎣⎦ Fig 3 Signals ( )() 1xt and 2xt Parameters A, T and ω 0 are chosen in the program „lfft” by the user , (Xω and(Xω Compute and plot the CTFT ( ))) 1Xω13 Compute the DFT of the discrete-time signals that result by sampling the three signals, with 8, 16 and 32 samples per period T Use the Matlab program „lfft” Compare the last two members of eq (28) (see also (24)) for each sample Compute the DFT using 256 samples Compare the plots obtained using the program with the computed CTFT and find the maximum error 5 2 Compute the convolution of the signals from 5 1: () , (xtxt∗, (xtxt∗ ())())() 11xtxt∗1222 Plot the results Use 512 samples and compute their convolution using the program “lfft” Compare the plots obtained by computation and the ones from Matlab 5